Spin-dependent localized Hartree-Fock density-functional approach for the accurate 
treatment of inner-shell excitation of close-shell atoms 
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Abstract 

We present a spin-dependent localized Hartree-Fock (SLHF) density-functional approach for the treatment of the inner-shell 
excited-state calculation of atomic systems. In this approach, the electron spin-orbitals in an electronic configuration are 
obtained first by solving Kohn-Sham (KS) equation with SLHF exchange potential. Then a single-Slater-determinant energy 
of the electronic configuration is calculated by using these electron spin-orbitals. Finally, a multiplet energy of an inner-shell 
excited state is evaluated from the single-Slater-determinant energies of the electronic configurations involved in terms of 
Slater's diagonal sum rule. This procedure has been used to calculate the total and excitation energies of inner-shell excited 
states of close-shell atomic systems: Be, B^, Ne, and Mg. The correlation effect is taken into account by incorporating the 
correlation potentials and energy functionals of Perdew and Wang's (PW) or Lee, Yang, and Parr's (LYP) into calculation. 
The calculated results with the PW and LYP energy functionals are in overall good agreement with each other and also 
with available experimental and other ab initio theoretical data. In addition, we present some new results for highly excited 
inner-shell states. 

PACS numbers: 31.15.Ew, 32.80.Wr, 32.80.Rm 
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I. INTRODUCTION 

Density functional theory (DFT) has been widely 
applied to many areas in theoretical physics and chem- 
istry as a powerful ab initio approach for the calculation 
of ground-state properties of many-electron systems due 
to its computational simplicity and efficiency 0, Q . The 
basic equation of DFT is Kohn-Sham (KS) equation Pj 
and the key part in KS equation is exchange-correlation 
(XC) potential 0. 

DFT with a traditional XC potential obtained from 
uniform electron gas, such as local density approxima- 
tion (LDA) IdNland generalized gradient approximation 
(GGA) ||J,UMl3, is a ground-state approach. Because 
of incomplete cancellation of spurious self-interactions in 
the conventional DFT using LDA or GGA [1 0, B 
and the inherent degeneracy (due to the use of spin 
and angular-momentum independent local potentials), 
the differences of the KS energy eigenvalues of unoc- 
cupied and occupied orbitals are not rigorously defined 
as excitation energies. However, the KS energy eigen- 
values can serve as good zeroth-order excited-state en- 
ergies provided they are obtained by solving KS equa- 
tion with a high-quality XC potential A number 
of theoretical methods have been developed by adopt- 
ing this point of view In particular, density 
work- functional approach (WF) 15', Tsl IT^ lisj . open- 
shell localized Hartree-Fock (LHF) density-functional ap- 
proach 0, 0, 0, ll^ . and multireference LHF density- 
functional approach |20ll2ll|. etc., have been successfully 
used to calculate excited-state properties of atomic and 
molecular systems. 

Recently, an exchange (X)-only LHF density- 
functional theory has been proposed and successfully ap- 
plied to ground-state calculations of atomic and molecu- 
lar systems In this X-only DFT, the exchange po- 



tential in the KS equation is a LHF exchange potential 
derived under the assumption that X-only KS determi- 
nant is equal to the Hartree-Fock (HF) determinant. We 
have recently extended this approach to excited states 
of atomic and molecular systems by assuming that the 
X-only KS determinant is also equal to the HF determi- 
nant for excited states |2^. Based on this postulate we 
have developed a spin-dependent localized Hartree-Fock 
(SLHF) density-functional approach for excited-state cal- 
culation of atomic and molecular systems 22]. In this 
approach, the exchange potential in the KS equation is 
an exact nonvariational SLHF exchange potential con- 
structed for both the ground and excited states. The 
SLHF potential is an analogue of the LHF potential. It 
is self-interaction free and exhibits the correct long-range 
behavior. Further, the SLHF potential requires the use 
of only the occupied orbitals and is dependent of the or- 
bital symmetry of the state. This approach associating 
with Slater's diagonal sum rule [23|| has been successfully 
used to calculate singly, doubly, and triply excited states 
of valence electrons of He- and Li-like ions '2^1 with ac- 
curate results. 

In this paper, we extend the SLHF density-functional 
approach to inner-shell excited states of atomic systems. 
We compute the total and excitation energies of inner- 
shell excited states of close-shell atomic systems: Be, B"*", 
Ne and Mg. In the calculation, the correlation potentials 
and energy functionals proposed by Perdew and Wang 
(PW) and by Lee, Yang, and Parr (LYP) are used 
to take into account the electron correlation effect. We 
will show that the calculated results are in overall good 
agreement with available theoretical and experimental 
data, demonstrating that the SLHF density-functional 
approach can provide a simple and computationally effi- 
cient approach for the accurate calculation of inner-shell 
excited states of close-shell atomic systems within DFT. 
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Finally, we also present some new results for the highly 
excited inner-shell states for the first time. 



II. THEORETICAL METHOD 

The SLHF density-functional approach has been dis- 
cussed in Ref. in detail and is outlined in this section 
for convenience. 

In spin-dependent density-functional approach, a spin- 
orbital ifia- (r) of the ith electron with spin a {a = a and 
f3 for spin-up and spin-down, respectively) and its orbital 
energy Ei^ are determined by the KS equation 



is the Slater potential |23| and 



where. 



Hair)ipta (r) = Eiaipio- (r) , 



(1) 



(2) 



is the KS Hamiltonian and 

Vf (r) = Ve^t (r) + Vh (r) + 14,, (r) , (3) 

is the local effective potential. In Eq. ((SJ, Vext (r) is the 
external potential, Vh (r) is Hartree potential (classical 
Coulomb electrostatic potential between electrons), and 
Vxca (r) is the XC potential. 

For a given atomic system, the external potential 
Vext (r) is known exactly. The Hartree potential Vh (r) 
is given by 



(4) 



where, p (r) = pa (r) + pp (r) is the total electron den- 
sity and Pa- (r) (for a ~ a and /?) is the spin-dependent 
electron density defined by 



P<T (r) ^^w,a\^^a (r)|^ 



(5) 



Here is the number of electrons with spin cr and Wi^ 
is the occupied number of electrons in the spin-orbital 
Via (r)- 

The XC potential can be decomposed into the ex- 
change potential Vxa (r) and the correlation potential 
Vca (r)- In the SLHF density-functional approach, the ex- 
change potential is a SLHF exchange potential V^^^^ (r) . 
It is given by 



Vl^^^{r) = Vl{r) + V^^Ar), 



where. 



Vxlir) 



1 



(6) 



(8) 



is a correction to Slater potential. In Eqs. (O and ^ 
j'^j (r) and Q^- are defined by 



lij (r) = ip.,air)ipja{r), 



and 



SLHF 



Vxa Vi> 



(9) 



(10) 



where, V^J" is a nonlocal exchange operator of the form 
of HF exchange potential but constructed from KS spin- 
orbitals. 

The SLHF exchange potential determined by Eqs. 
(|10|l has two arbitrary additive constants. The physical 
orbitals can only be obtained by the use of appropriate 
constants in the exchange potential To settle down 
the constants so as to pick up the physical orbitals, it 
is required that the highest-occupied-orbital No- of each 
spin a does not contribute to the correction term V^^(r). 
In this case, the correction term V^^i''^) decays exponen- 
tially, the SLHF exchange potential behaves asymptoti- 
cally as Slater potential and thus approaches to — 1/r at 
long range ,10j. 

In atomic systems, an electron spin-orbital is charac- 
terized by three quantum numbers n, I, and cr, where 
n and I are the principal quantum number and orbital 
angular momentum quantum number of the electron, re- 
spectively. In the spherical coordinates, the spin-orbital 
(fiia- (r) of an electron with quantum numbers n, I, and a 
can be expressed by 



ifiia (r) = Yir^ 



(11) 



where, Rnia-{r) is the radial spin-orbital, Yim{0,(l)) is the 
spherical harmonic, m is the azimuthal quantum number, 
and I is a set of quantum numbers apart from spin a of the 
spin-orbital. The radial spin-orbital Rnia-{r) is governed 
by radial KS equation. 



2dr^ 



1(1 + I) 
2r2 



„cff 



(r) 



where v^{r) is the radial effective potential given by 
vf{r) = Vext (r) + VH (r) + t;SLHF (^) ^ (^3) 

In Eq. (Unjl, Vext {r), vh (r), wf^"^ (r), and Vca {r) are 
the radial external potential, radial Hartree potential, 
radial SLHF exchange potential, and radial correlation 
potential, respectively. 

For an atomic system with a nuclear charge Z, the ex- 
ternal potential is the Coulomb potential between elec- 
tron and nucleus 



Vext (r) 



Z 



(14) 
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In central-field approach, the radial Hartree potential 
is calculated from 



vh (r) ^ An I —p{r')r''^dr', 



(15) 



where, r> is the larger of r and r', p{r) ~ Pair) + pf3{r) is 
the spherically averaged total electron density, and po- (r) 
(cr = a or /3) is the spherically averaged spin-dependent 
electron density given by 



1 

nl 



Rn 



(16) 



Here the symbol v„ stands for a set of quantum num- 
bers for summation and the sum is performed over all 
the occupied spin-orbitals with spin a. This expression 
is accurate for spherically symmetric (close-shell) states, 
but it is only an approximation for non-spherically sym- 
metric (open-shell) states. It may induce an error when 
it is used to evaluate the energy of a non-spherically sym- 
metric state. However, the error is negligible compared 
to the order of calculated multiplet splitting 1 1 . 
The radial SLHF exchange potential is given by 



(r) = vl, (r) -f (r) , 



(17) 



where, 



fa Va 



via ir) 



■EE 



47rp<j(r) 



nlm n' I'm' 

is the radial Slater potential and 



''^xa- (r) 47^p^(^^-) E/ E/ '^nlm,n'l'va> {''')■ 



(19) 



nlm n'l'm' 



is a correction to the radial Slater potential. The matrix 
elements sli^^^,{r) and (r) in Eq. ^ are 

given in Ref. |22| . 

To calculate electron spin-orbital, the Legendre gen- 
eralized pseudospectral (LGPS) method is used to 
discretize the radial KS equation ifT^ . This method asso- 
ciated with an appropriate mapping technique can over- 
come difficulties due to singularity at r = and long-tail 
at large r of Coulomb interaction and thus provides a 
very effective and efficient numerical algorithm for high- 
precision solution of KS equation. Using the electron 
spin-orbitals of an electronic configuration, a single Slater 
determinant for a specific electronic state is constructed 
and its total energy calculated. The total energy is a 
sum of non-interacting kinetic-energy Ek, external-field 
energy E^xt, Hartree energy Eh, exchange energy E^, 
and correlation energy Ec- The values of Ek, E^xt, Eh, 
and Ex are evaluated by 



Ek 



/3 i'CT /. 

E E ^"'"^ / 

ry=a nl 



1(1 + 1) 
2r2 



Eext = 47r / Vext (r) p (r) r^dr, 



Vlrr. 



2 ^ 
n 



i' nla,i 



'I'cr'^ 



and 



Ex 



E 

n 



^In 



k 

nla.n'V ' 



.'5 a 



(21) 



(22) 



(23) 



where, H represents a collection of all the quantum num- 
bers involved, the matrix elements rjimi'm'^ -E!^ia n'Va' ^ 



and Gz 



are given in Ref. |22 



For a multiplet state that can be described completely 
by a single Slater determinant, the energy is calculated 
directly from the single Slater determinant. For a multi- 
plet state that cannot be represented by a single determi- 
nant, the energy can be calculated by means of Slater's 
diagonal sum rule [2^ . According to this rule, a sum over 
single-Slater-determinant energy E{Y}i) of determinant 
Di from an electron configuration equals to a weighted 
sum over multiplet energy £'(Mj) of multiplet state Mj 
involved in the same electron configuration, namely. 



Y,Em^Y.^jE{M,), 



(24) 



Rnla (r) 

(20) 



where, the weight dj is the times that the multiplet state 
Mj appears in all the single Slater determinants. Similar 
procedures have been employed in recent excited-state 
calculations [11 El 111 . 



III. RESULTS AND DISCUSSION 

The procedure described in the preceding section is ex- 
tended to calculate the total energies (E) and excitation 
energies (Ai?) of inner-shell excited states of closed-shell 
atomic systems: Be, B+, Ne, and Mg. In the calcu- 
lations, the correlation effect, which is characterized by 
the correlation potential Vca (f) and correlation energy 
Ec, is taken into account through the correlation poten- 
tials and energy functional of Perdew and Wang (PW) 
and of Lee, Yang, and Parr (LYP) Q, respectively. 
The results obtained with these two correlation energy 
functionals are listed in columns PW and LYP in the 
following tables, respectively. For simplicity, hencefor- 
ward we use an abbreviation (nl) ^ {n'V) to represent an 
electronic configuration for an electron in inner-shell {nl) 
being excited to sub-shell [n'l'), unless otherwise speci- 
(^d. For instance, an abbreviation ls~^2p represents an 
electronic configuration of inner-shell excitation ls2s^2p 
of Be. 
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A. Be 

In TABLE we present the total energies and exci- 
tation energies from our calculations for inner-shell ex- 
cited states ls~^np ^'^P (n = 2 ~ 8) of Be. For com- 
parison we also list in this table the theoretical results 
of density work- functional approach (WF) J2^, saddle- 
point complex-rotation approximation (SPCR) [27i l28l|. 
and R-matrix method in the close-coupling approxi- 
mation (RMCC) ||29| and experimental results (Exp.) 
p50, 31, 32, 33]. For the total energies, the maximum rel- 
ative discrepancies of our PW and LYP results are 0.10% 
and 0.25% to the WF results, and 0.27% and 0.39% to the 
SPCR results. For excitation energies, the maximum rel- 
ative deviations of our PW and LYP results to the exper- 
imental results are 0.37% and 0.90%, respectively, while 
the maximum relative discrepancies of the WF, SPCR, 
and RMCC results to the experimental results are 0.41%, 
0.05%, and 0.14%, respectively. This demonstrates that 
both the total energies and excitation energies from our 
calculations agree well with the experimental and other 
theoretical results. On the other hands, due to overesti- 
mation of LYP energy functional to correlation energies 
of atomic systems with smaller Z |^ the LYP re- 
sults are a little bit worse than the PW results. It is also 
shown that our PW results for excitation energy are a 
little bit better than the WF results. 

B. B+ 

To explore the feasibility of the approach to inner-shell 
excitation of atomic ions, we also apply the procedure to 
inner-shell excited-state calculation of B"*". The total en- 
ergies and excitation energies of inner-shell excited states 
Is'^np ^'^P (n = 2 ~ 8) are given in TABLE Hll along 
with theoretical results of Dirac-Fock method (DF) |3J| 
and available experimental results [13 ■ For excitation 
energies, the relative deviations of our PW and LYP re- 
sults are less than 0.32% and 0.64% to the DF results, 
and less than 0.29% and 0.60% to the experimental re- 
sults. This demonstrates again that our results are in 
good agreement with both the experimental and other 
theoretical results, and the PW results are a little bit 
more accurate than the LYP results for this atomic ion 
having smaller Z. 



results of density work- functional ap pro ach (WF) [2^ . 
configuration- interaction model (CI) [33, and Hartree- 
Fock method (HF) l3fl , a nd experimental results (Exp.) 

For the total energies of excited states ls~^ns ^'^S 
(n = 3 ~ 8) given in TABLE IIIII the relative deviations 
of our PW and LYP results to the WF results are not 
more than 0.03% and 0.02%, respectively. This demon- 
strates that the PW energy functional has almost the 
same precision as the LYP energy functional in calcula- 
tion of total energies of these inner-shell excited states. 
For excitation energies, the relative discrepancies of our 
PW and LYP results to the experimental result are less 
than 0.13% and 0.02%, respectively. This indicates that 
the LYP results for the excitation energy are better than 
the PW results. Since the total energies from calcula- 
tion with PW energy functional are very close to those 
with LYP energy functional, the larger discrepancy of the 
PW results for excitation energy mainly comes from the 
ground-state energy, which is -128.8952 (a.u.) from the 
calculation with PW energy functional. This value is dif- 
ferent from that with LYP energy functional -128.9331 
(a.u .) and that obtained from Hartree-Fock energy 
|45| plus correlation energy ||46i] -128.937 (a.u.). In ad- 
dition, the maximum relative discrepancies of the exci- 
tation energies from WF calculation and CI calculation 
are 0.08% and 0.03% to the experimental results, respec- 
tively. This illustrates that our LYP results are very close 
to the CI results and better than the WF results. 

For the total energies of inner-shell excited states 
Is-^np ^-^P (n = 3 - 8) given in TABLE |IV| the rel- 
ative deviations of our PW and LYP results to the WF 
results are not more than 0.02% and 0.02%, respectively. 
This implies that the PW energy functional has the same 
precision as the LYP energy functional in the total energy 
calculation of these states. For the excitation energies, 
the maximum relative discrepancies of our PW and LYP 
results to the experimental results are 0.14% and 0.02%, 
while the maximum relative deviations of the WF and CI 
results to the experimental results are 0.07% and 0.02%, 
respectively. This demonstrates that the LYP results for 
excitation energy are a little bit more accurate than the 
PW results for these inner-shell excited states. It is also 
shown that the LYP results are again very close to the 
CI results and a little bit better than the WF results. 



C. Ne 

We present in TABLE UTIl and TABLE |IV| the total 
energies and excitation energies from our calculations 
for inner-shell excited states ls~^ns ^'^S (n = 3 ~ 8) 
and Is^^np ^''^P (ti = 3 8) of Ne, respectively. We 
also present in TABLE IIVI the total energy of inner- 
shell excited state ls2s^2p^ '^S of Ne+ and ionization 
energy of an inner-shell Is electron of Ne. For com- 
parison we also show in these tables the theoretical 



D. Mg 

For Mg, we have computed the total energies and 
excitation energies for inner-shell excited states 2p^^ns 
i^3p (n = 4 ~ 8), 2s-^np ^-^P (n = 3 ~ 8), 2s-^ns ^-^S 
(n = 4 - 8), Is-^np ^-^P (?i = 3 - 8), and Is-^ns ^-^S 
(„ = 4 - 8). The results are shown in TABLE |Vl to TA- 
BLE IIXI respectively. For comparison, the theoretical re- 
sults from configuration-interaction calculation with im- 
proved and optimized orbitals (CIIOO) j43| and exper- 
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TABLE I: Total energies (E) and excitation energies {AE) of inner-shell excited states ls~^np ^'^P (n = 2 ~ 8) of Be. The 
ground state energies obtained from calculation with PW and LYP correlation potentials and energy functionals are —14.6575 
(a.u.) and —14.6686 (a.u.), respectively. Here 1 a.u.=27.2116 eV is used. 
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"PW results, ''LYP results, =|2^, ''I^IH, f^^, »|3l, and ''I3I. 



TABLE II: Total energies (E) and excitation energies (AE) of inner-shell excited states ls~^np ^'^P (n = 2 ~ 8) of B+. The 
ground state energies obtained from calculation -with PW and LYP correlation potentials and energy functionals are —24.3284 
(a.u.) and —24.3432 (a.u.), respectively. 
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217.1886 


217.8745 






ls-^8p 
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16.3362 


16.3256 


217.4803 


218.1725 








IP 


16.3359 


16.3254 


217.4895 


218.1767 







"P"W results, ""LYP results, =[3, and "^Isl. 



TABLE III: Total energies (E) and excitation energies (AE) of inner-shell excited states ls~^n3 ^'^S (n = 3 ~ 8) of Ne. The 
ground state energies from calculations "with PW and LYP correlation potentials and energy functionals are -128.8952 (a.u.) 
and -128.9331 (a.u.), respectively. 
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97.1495 
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864.6583 


865.37 


865.1 865.32 
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97.0335 


97.0348 




867.0081 


868.0038 
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97.0319 
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867.0516 


868.0636 




868.21 
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96.9999 


97.0008 




867.9208 


868.9290 








Is 


96.9993 


97.0000 




867.9382 


868.9507 




869.06 




ls-^6s ^S 


96.9853 


96.9863 




868.3189 


869.3235 








^S 


96.9850 


96.9859 




868.3276 


869.3344 




869.44 






96.9776 


96.9788 




868.5279 


869.5276 








Is 


96.9774 


96.9786 




868.5328 


869.5331 








Is-^Ss ^S 


96.9731 


96.9744 




868.6522 


869.6474 








^S 


96.9729 


96.9743 




868.6555 


869.6501 









"P"W results, '"LYP results, =|^|3l, and -^Issl. 
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TABLE IV: Total energies (E) and excitation energies (A_E) of inner-shell excited states Is ^np ^'^P (n = 3 ~ 8) of Ne and 
ls2s^2p^ of Ne+. 



States 




-B 




(a.u.) 




AB 








(cV) 




Present 


work 


Other theory 


Present 


work 


Other 


theory 




Exp. 




PW° 


LYF" 


WF'= 


PW" 


LYP" 


WF= 


CI" 


HB" 


W' 


ADC 




■ip 


97.0762 


97.0766 


97.0982 


865.8456 


866.8663 


866.4244 














97.0743 


97.0736 


97.0950 


865.8984 


866.9480 


866.5115 


867.18 


867.05 


867.13 


867.12 




3p 


97.0147 


97.0151 


97.0330 


867.5194 


868.5398 


868.1986 












ip 


97.0140 


97.0141 


97.0318 


867.5377 


868.5671 


868.2312 


868.70 


868.68 


868.77 


868.69 


Is^'-bp 


3p 


96.9914 


96.9919 


97.0098 


868.1540 


869.1712 


868.8299 












IP 


96.9910 


96.9915 


97.0095 


868.1627 


869.1820 


868.8462 


869.32 


869.23 


869.37 


869.27 


Is^^Gp 


3p 


96.9806 


96.9814 


96.9900 


868.4457 


869.4569 


869.1238 












IP 


96.9805 


96.9812 


96.9988 


868.4506 


869.4623 


869.1347 


869.58 


869.63 


869.65 


869.56 


Is^^Tp 


3p 


96.9747 


96.9756 




868.6071 


869.6147 














IP 


96.9746 


96.9755 




868.6092 


869.6174 




869.79 






869.73 


Is^^Sp 


3p 


96.9710 


96.9721 




868.7088 


869.7099 














IP 


96.9710 


96.9720 




868.7088 


869.7127 




869.87 








ls2s^2p'' 




96.9612 


96.9629 




868.9741 


869.9603 


869.6898 


870.15 


870.10 


870.17 


870.17 



"PW results, 'LYP results, '=|2^, (see also =|33|3l, -^^3. and "[S^ (see also ISSllIol V 



imental results (Exp.) [isl l49l | are also shown in these 
tables. 

For inner-shell excitation of Mg, the excited states in- 
volving in excitation of a 2p electron, 2p~^ns ^'^P, are the 
only ones that experimental excitation energies are avail- 
able, as shown in TABLE0 For these states, the excita- 
tion energies from our calculation are in good agreement 
with the experimental results. The relative deviations of 
our PW and LYP results to the experimental results are 
not more than 0.38% and 0.36%, respectively. Thus the 
PW energy functional has the same precision as the LYP 
energy functional in the excitation energy calculation of 
these inner-shell excited states. Apart from an excited 
state 2p~^4s our results also agree well with the CI- 
100 results. The maximum relative discrepancies of our 
PW and LYP results to the CIIOO results are 0.58% and 
0.60%, respectively. The larger discrepancies of our re- 
sults to the CIIOO results for the state 2p~^4s ^P are 
caused by the fact that the excitation energy from the 
CIIOO calculation is too small. It is shown from TA- 
BLE that the CIIOO result is much smaller than the 
experimental data while our excitation energies are very 
close to the experimental results. 

For inner-shell excited states relevant to excitation of a 
2s electron, the excitation energies from our calculations 
with both PW and LYP energy functionals are larger 
than the CIIOO results, as shown in TABLEEHand TA- 
BLE IVIll For the excited states 2s~^np '^P, except an ex- 
cited state 2s~^3p "^P, the relative deviations of our PW 
and LYP results to the CIIOO results are less than 0.70% 
and 0.62%, respectively. For the excited state 2s~^3p '^P 
our PW and LYP results are much larger than the CI- 
IOO result. For the excited states 2s~^ns ^S, apart from 
an excited state 2s~^4s ^S, the relative discrepancies of 
our PW and LYP results to the CIIOO results are not 
more than 0.41% and 0.25%, respectively. For the ex- 
cited state 2s~^4s our PW and LYP results are again 
much larger than the CIIOO results. 

For excitation energies of excited states ls~^np ^P and 



Is-^ns given in TABLE IVlTTl and TABLE El our 
PW and LYP results are smaller and larger than the 
CIIOO results, respectively. For excited states ls~^np 
^P the relative deviations of our PW and LYP results 
to the CIIOO results are less than 0.02% and 0.08%, 
respectively. For excited states \s~^ns ^S, apart from 
a state ls~^4s ^S', the relative discrepancies of our PW 
and LYP results to the CIIOO results are not more than 
0.04% and 0.03%, respectively. For the state Is'^As 
our LYP result is larger than the CIIOO result by 0.11%. 



IV. CONCLUSIONS 

In summary, the procedure we have developed for 
excited-state calculation based on SLHF density func- 
tional approach and Slater's diagonal sum rule has been 
extended to the treatment of inner-shell excited states of 
atomic systems. In this procedure, electron spin-orbitals 
in an electronic configuration are obtained first by solving 
the KS equation with the exact SLHF exchange potential. 
Then a single-Slater-determinant energy of the electronic 
configuration is calculated by using these electron spin- 
orbitals. Finally, a multiplet energy of an excited state is 
evaluated from the single-Slater-determinant energies of 
the electronic configurations involved in terms of Slater's 
diagonal sum rule. In this procedure, the key part is the 
SLHF exchange potential. This potential qualifies for 
inner-shell excited-state calculation because it provides 
a potential with free self-interaction, correct long-range 
behavior, and symmetry dependence of atomic state. We 
have applied this procedure to the calculations of to- 
tal energies and excitation energies of inner-shell excited 
states of close-shell atomic systems: Be, B+, Ne, and 
Mg. In the calculation, the generalized pseudospectral 
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TABLE V: Total energies (E) and excitation energies (AE) of inner-shell excited states 2p~^ns ^'^P (n = 4 ~ 8) of Mg. The 
ground state energies from calculations with PW and LYP correlation potentials and energy functionals are -200.0204 (a.u.) 
and -200.0744 (a.u.), respectively. 



States 


-E 


(a.u.) 




AE 




(eV) 




Present 


work 


Present 


work 


Other theory 


Exp. 




PW 


LYP" 


PW 


LYP" 


CIIOO'= 


MZ" 


NE" 


2p-Ms 


3p 


198.0039 


198.0557 


54.8714 


54.9323 


53.72 


54.801 


54.801 




IP 


198.0022 


198.0538 


54.9190 


54.9849 




55.065 


55.065 




3p 


197.9509 


198.0044 


56.3133 


56.3272 


55.99 


56.278 


56.280 






197.9504 


198.0038 


56.3294 


56.3457 




56.544 


56.545 


2p"i6s 


3p 


197.9313 


197.9852 


56.8469 


56.8499 


56.56 


56.777 


56.785 




IP 


197.9312 


197.9849 


56.8505 


56.8586 




57.039 




2p"^7s 


3p 


197.9217 


197.9759 


57.1095 


57.1038 










IP 


197.9215 


197.9757 


57.1142 


57.1087 




57.302 


57.305 


2p"i85- 


3p 


197.9163 


197.9707 


57.2573 


57.2464 










IP 


197.9161 


197.9706 


57.2603 


57.2478 






57.456 



"PW results, ""LYP results, and °l4l. 



TABLE VL Total energies (E) and excitation energies (A_E) TABLE VIII: Total energies (E) and excitation energies (A_E) 
of inner-shell excited states 2s~^np ^'^P (n = 3 ~ 8) of Mg. of inner-shell excited states ls~^np ^'^P (n = 3 ~ 8) of Mg. 





-E 


(a.u.) 


AE 




(eV) 




-E 


(a.u.) 


AE 




(eV) 


States 


Present 


work 


Present 


work 


Other theory 


States 


Present 


work 


Present 


work 


Other theory 




PW" 


LYP" 


PW" 


LYP° 


CIIOO'^ 




PW" 


LYP" 


PW" 


LYP" 


CIIOO" 


2s"^3p ^P 


196.6182 


196.6721 


92.5785 


92.5820 


91.72 


ls~^3p ■'P 


152.1281 


152.1433 


1303.2275 


1304.2822 


1303.25 


ip 


196.6101 


196.6630 


92.8005 


92.8288 




ip 


152.1207 


152.1353 


1303.4275 


1304.5002 




2s"^4p ^P 


196.4612 


196.5181 


96.8526 


96.7718 


96.18 


Is^Up ^P 


151.9671 


151.9851 


1307.6080 


1308.5868 


1307.86 


ip 


196.460 


196.5170 


96.8820 


96.8026 




ip 


151.9659 


151.9839 


1307.6412 


1308.6192 




2s"^5p ^P 


196.4279 


196.4871 


97.7574 


97.6164 


97.11 


ls"^5p ^P 


151.9334 


151.9536 


1308.5256 


1309.4453 


1308.80 


ip 


196.4275 


196.4865 


97.7677 


97.6317 




ip 


151.9329 


151.9533 


1308.5375 


1309.4521 






196.4138 


196.4738 


98.1403 


97.9775 




Is^'^Gp ^P 


151.9191 


151.9401 


1308.9131 


1309.8108 




ip 


196.4136 


196.4735 


98.1463 


97.9863 




ip 


151.9189 


151.9399 


1308.9196 


1309.8179 




2s"i7p "P 


196.4065 


196.4669 


98.3403 


98.1667 




Is'^lp ^P 


151.9117 


151.9331 


1309.1150 


1310.0018 




ip 


196.4064 


196.4667 


98.3441 


98.1721 




ip 


151.9116 


151.9329 


1309.1188 


1310.0062 




2s"^8p ^P 


196.4022 


196.4628 


98.4581 


98.2782 




ls"^8p ^P 


151.9073 


151.9290 


1309.2344 


1310.1145 




ip 


196.4021 


196.4626 


98.4603 


98.2820 




ip 


151.9072 


151.9289 


1309.2366 


1310.1172 





"PW results, ''LYP results, and "PW results, ''LYP results, and "Wi. 



TABLE VII: Total energies {E) and excitation energies (AS) TABLE IX: Total energies {E) and excitation energies (Ai5) 
of inner-shell excited states 2s~^ns ^'^ S (n = 4 ~ 8) of Mg. of inner-shell excited states ls~^ns ^'^ S (n = 4 ~ 8) of Mg. 



States 


~E 


(a.u.) 


AE 




(eV) 


States 


-E 


(a.u.) 


AE 




(eV) 


Present 


work 


Present 


work 


Other theory 


Present 


work 


Present 


work 


Other theory 


PW" 


LYP" 


PW" 


LYP" 


CIIOO'^ 


PW" 


LYP" 


PW" 


LYP" 


CIIOO'' 


2s"Ms 


196.4919 


196.5499 


96.0167 


95.9065 


94.65 


ls"-'4s '■'S 


151.9980 


152.0171 


1306.7666 


1307.7160 


1306.29 




196.4873 


196.5447 


96.1419 


96.0502 






151.9957 


152.0144 


1306.8281 


1307.7895 




2s"^5s 


196.4387 


196.4984 


97.4635 


97.3079 


97.07 


ls"^5s 


151.9442 


151.9649 


1308.2314 


1309.1362 


1308.77 




196.4372 


196.4966 


97.5060 


97.3585 






151.9434 


151.9640 


1308.2521 


1309.1618 






196.4191 


196.4792 


97.9985 


97.8325 


97.64 




151.9243 


151.9454 


1308.7716 


1309.6660 


1309.35 




196.4183 


196.4783 


98.0181 


97.8565 






151.9240 


151.9450 


1308.7808 


1309.6780 




2s"^7s 


196.4094 


196.4698 


98.2616 


98.0872 






151.9146 


151.9360 


1309.0366 


1309.9226 






196.4090 


196.4693 


98.2720 


98.1003 






151.9144 


151.9358 


1309.0421 


1309.9297 




2s-'^Ss 


196.4039 


196.4645 


98.4105 


98.2301 






151.9091 


151.9307 


1309.1868 


1310.0663 




Is 


196.4037 


196.4642 


98.4176 


98.2382 




Is 


151.9089 


151.9306 


1309.1901 


1310.0706 





"PW results, ''LYP results, and '=^. "PW results, ''LYP results, and "iM- 



method with nonuniform grids is used for optimal dis- 
cretization of the spatial coordinates, allowing accurate 
and efficient treatment of the KS equation and the spin- 
orbital energies for both the ground and excited states. 
The correlation effect is taken care of by incorporating 
the PW and LYP correlation potentials and energy func- 
tionals into calculation. The results from our calculations 
with LYP and PW energy functionals are in overall good 
agreement with each other and also with the available 



more sophisticated ah initio theoretical results and exper- 
imental data. The maximum relative discrepancy of our 
calculated excitation energies to the available experimen- 
tal results is not more than 0.90%, demonstrating that 
the SLHF density-functional approach is capable of pro- 
viding a powerful and computationally efficient scheme 
for accurate inner-shell excited-state calculation of close- 
shell atomic systems within DFT. Extension of the SLHF 
density-functional approach to open-shell atomic systems 
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is in progress. 
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